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Abstract 

In this paper, we show how to improve and extend the Integration by Fractional Expansion technique 
(IBFE) by applying it to certain families of scalar massive Feynman diagrams. The strategy is based on 
combining this method together with the Integration by Parts technique (IBP). In particular, we want to 
calculate certain Feynman diagrams which have a triangle loop as a subgraph. The main idea is to use IBP 
in this subgraph in order to simplify the topology of the original diagram in which it is immersed, using 
then, in a second step, the IBFE technique. The result we have obtained, after the application of both 
techniques, represents a simplification in the complexity of the solution, compared with having used only the 
IBFE technique. 
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1 Introduction 



Since the beginning of the development of quantum field theory, the evaluation of Feynman diagrams has 
been one of the most important theoretical challenges. The development of methods leading to the analytical or 
numerical evaluation of loop integrals represents an important problem both from the point of view of physics, 
as well as from the perspective of advanced mathematical methods. 

Concerning physics, this research is relevant for perturbative calculations of higher orders, needed to test our 
present theoretical models when their predictions are compared to experimental results from present and future 
colliders and/or from other high energy physics experiments. The amount of experimental information at present 
is so huge, and with such a high precision, that the numerical calculations are almost not enough to achieve 
the needed accuracy in order to compare to the experimental results. Therefore, it is extremely important to 
be able to optimize the present existing integration techniques and/or to generate new mathematical techniques 
that may allow us to face up the calculation of Feynman diagrams with a higher number of loops, or that may 
involve also the presence of multiple scales. 

In this perspective, we consider the Integration by Fractional Expansion Technique (IBFE) as a method that 
allows the calculation of certain families of Feynman diagrams in a simple and direct way, getting solutions for 
the general case, i.e. arbitrary indices for propagators, in terms of generalized hypergeometric functions. By 
indices we understand here, as usual, the power degree of the propagators. In some cases, for some special 
values of these indices, it is possible to carry out the sum of these series in terms of conventional functions which 
simplifies the subsequent analysis of the final solution. 

The IBFE technique, in general aspects, is similar to the integration technique based on the Mellin-Barnes 
representation of the diagram [TJ |2] . Other authors have called this technique NDIM (Negative Dimensional 
Integration Method) [3j 0J [5] , which is based on the work by Halliday and Ricotta [6] . These authors make an 
analytical continuation of the D dimension, working in the frame of dimensional regularization (D — 4 — 2e) 
into negative values of D, which justifies the name of the method. However, we prefer to use the name IBFE 
due to the reasons explained in [7]. 

It is possible to extend the range of validity of IBFE by considering a new kind of diagrams that may be 
evaluated by means of this technique. In particular, in this article we will show that with the help of the 
Integration by Parts (IBP) method [5], it turns out that certain families of diagrams have more simple solutions 
compared with those that appear when only the IBFE method is used. It is important to remark that the 
existence of these simple solutions allows us to increase the degree of difficulty of the Feynman diagrams that 
can be handled with the combined IBFE + IBP techniques. We could also be able to consider diagrams with a 
certain particular mass distribution in the propagators. 

To justify this methodology, we will concentrate on a certain class of Feynman diagrams that involve in their 
geometry a one-loop triangle subgraph, and such that we may introduce the IBP technique through the so-called 
triangle identity. The main idea is to use IBP to transform the geometry of the diagram. The elimination of 
certain propagators of the original diagram may produce now a diagram which will be ideal for the application 
of the IBFE method, as was pointed out in [7j. Basically, it turns out that triangles are transformed into bubbles 
under some conditions determined by the indices of the propagators. 

This work has been organized in the following way: In Sec. II, we remind the formula associated to the 
triangle identity. This rule emerges when the IBP technique is applied to a one-loop diagram with three external 
lines. Here we have assumed a scalar massless theory. 

In Sec. Ill, and in order to justify the combined use of IBFE + IBP techniques, we will evaluate a simple 
diagram, the one-loop triangle diagram with one massive propagator with mass to. To simplify the problem even 
more , we have assumed that two of the external lines are on mass shell with to = 0. In this way, the diagram 
will be characterized by only two invariants {to 2 , s}. The above mentioned simplifications, as we will see later, 
do not imply a loss with respect to the general statements which are the goal of this section. 

The validation of the combined use of IBFE + IBP will be realized through a calculation of this diagram, 
using first only the IBFE technique. Later, the same diagram is evaluated with the combined techniques IBFE 
+ IBP. We will show that both techniques provide equivalent solutions, both for numerical comparisons for some 
specific values of the invariants of the problem, as well as for analytical comparisons in certain special limits, 
(to = or s = 0). 
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In Sec. IV, we present a more complex application. The idea is to specify the range of validity of the combined 
use of the IBFE + IBP techniques. In particular, we consider the topology associated to the radiative corrections 
to the propagator with two loops and five internal lines. This diagram can be easily evaluated by means of IBP 
when all of the propagators are massless with index one. In our case, we will extend the discussion by considering 
now this diagram with one massive propagator and allowing some arbitrary indices. The analytic solutions for 
this problem will be given in terms of generalized hypergeometric functions of the form q F q _i, whose argument 
(M 2 \ ( p 2 \ 

can be — =- or — =• , according to the relevant kinematical region; p denotes the momentum that flows 
V p 2 J \M 2 J 

into the diagram and which fixes the scale, and M is the mass in the massive propagator. The general solutions 
of Feynman diagrams will correspond always to hypergeometric series, i.e. to natural expansions around a zero 
value of the argument and with a convergence radius r = 1. The mathematical formalism, the relevant formulas, 
and applications are presented and discussed in [9[ [TOl QU HI] )• 



2 IBP and the triangle rule 

The well known triangle rule is obtained directly by applying the IBP method to the massless one-loop 
diagram with three external lines (how his rule is derived can be seen in the book by V. A. Smirnov |13|). The 
integral that represents the triangle diagram is the following: 

1 1 f d D q 1 1 1 m 

(p a ) 2ai {p P f a " J ^ D/2 1(p<* + <?) 2 P [(Pp - ?) 2 P (<zT 3 ' 1 ' 

where P a and Pp identify the independent momenta of the external lines. Then, associating each index 
{01,02,03} with his own propagator, we may introduce a graphical representation for the triangle identity 
once IBP has been applied to the previous integral. 




In this expression, D denotes the number of space dimensions in the integral. The indices ai (i = 1,2,3,4,5) 
correspond, as we said, to the powers of the propagators of the corresponding internal lines in the diagram. The 
indices 04 and a§ will be arbitrary if the corresponding lines are internal lines in a diagram that contains this 
triangle and their value will be zero when the lines are external. Obviously, the IBP technique goes beyond this 
identity. However, for our purpose this identity will be enough. 



3 Validation of the IBFE + IBP method 

In previous articles [3 [T5] , we have shown that the IBFE technique is a useful method for evaluating certain 
families of Feynman diagrams. In particular, we will use now this technique for handling the following test 
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topology: a one-loop triangle diagram where one of the propagators has a finite mass and the two external lines 
are massless. The result for this diagram will be compared with the result that emerges from the combined IBFE 
and IBP techniques. Explicitly, we will expand the massive propagator according to the IBFE prescription, then 
extracting out the mass from the diagram and getting finally a series that involves only non massive diagrams. 
Each of them may be reduced topologically, by means of the IBP approach, to a series of bubble diagrams. The 
IBFE technique will be applied then to each one of these bubble diagrams, obtaining the Multiregion Expansion 
(MRE) [7J of the complete diagram which allows us to obtain the solutions in terms of hypergeometric functions 
of the form q F q ^\. 

Once we get the solution of this diagram in both ways, we will make a analytical comparison in the limits 
m = and s = 0. 



3.1 IBFE solution 

The diagram to be evaluated is: 



G=Pn 




(3) 



such that the propagator associated to the index one has a mass m, whereas Pf = P$ = and P| = s [2J- For 
this diagram then, the equivalent configuration of invariants is described as follows: 




G = s 



We will start by writing the corresponding Schwinger parametrization for this topology which is given by 

X\X 2 



(4) 



G 



(-l)- p / 2 

nLi r ( a i) 



d!& exp (x\j 



exp 



Xi+X 2 + x 3 



(xi + x 2 + x 3 ) 



(5) 



where d~x* = dx\dx 2 dxz- Some algebra allows us to obtain the MRE [7] form of (JSJ): 



G 



(-ir D/2 , , ^«w, n2 n-=iA 



4>n u ...,n B (-to 2 )™ 1 (s) r 



n j= i r K) ni ,...,„ 5 

where the factor </> ni) ..,, n5 , defined in [7j is given as follows: 
The constraints are given by the identities 



r(f + n 2 ) 



1 



r(m + i)...r(»B + 1) 



(6) 



(7) 



Ai = + n 2 + n 3 + n 4 + n 5 ) 
A 2 = (ai +n l +n 2 + n 3 ) , 
A 3 = (a 2 +n 2 + n 4 ) , 
A 4 = (a 3 + n 5 ) . 



(8) 
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Each constraint, according to the IBFE technique [7] is associated to the replacement of the integral, after 
expanding the integrand, with the following symbol: 

dx x ai+ - +a2 - 1 = (ai + ... + a 2 ) (9) 
which also emerged as a result from multinomial expansions that obey the following rule [7] : 

(A, + ... + Aa) ±V = E - E ^.-...n- ^ + y + + (10) 

where the terms Ai (i = 1, a) and the exponent v are quantities that can take arbitrary values. 
3.1.1 Solution for arbitrary indices 

The properties of the general solutions for Fcynman diagrams imply the existence of two interesting kincmatical 
regions: 

Solution in the region 



< 1 



G 



where the factor r\ is defined as: 



ai + a 2 + a 3 



D 



«2 



D 

2 



| , »- ai -a 2 -a 3 r (oi + a 2 + a 3 - f ) r (f - a 2 - q 3 ) 

' 1 ; r(ai)r(f) 



(ii) 



(12) 



s 



< 1 



Solution in the region 



1 + ai + a 2 + a 3 - Z? , ai + a 2 + a 3 



1 + Oi + a 3 



_D 

2 



nr 

s 



(13) 



+ V X 2-Fl 



1 + a 2 



a 2 



1 + T - ai - a 3 



rn 
s 



where the different factors are: 



v = ( _ 1} -f (s) f - ai - 02 -a 3 r(a 1+ a 2 +a 3 -f)r(f - oi - a 3 ) r (f - a 2 - a 3 ) 



D 



2 =(-l)-T (sp 2 (-rn 2 ) 2 



3.1.2 Solution for unitary indices 



r (ai) T (a 2 ) T (L> - a : - a 2 - a 3 ) 
13 ai-a 3 r (ai + a 3 - f ) r - a 2 - a 3 ) 



r(oi)r(f -a 2 ) 



(14) 
(15) 



This situation is especially relevant, since when an arbitrary Feynman diagram is computed, the indices asso- 
ciated to the propagators are normally unitary, this case also being useful for numerical computations. Then, 
with D = 4 - 2e, the Eqs. (HI]) and (JT3]) are: 
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Solution in the region 



< 1 



c(^) = -(-ir (~m 2 ) 



_i_ e r(l + e) 



e(l-e) 



1F1 



1 + e , 1 

2 - e 



(16) 



Solution in the region — < 1 

s 



2e 



and the factors 



' 1 j [ ' e 2 r(l-2e) 



!Z =-(-ir( S )- 1 (-™ 2 ) 



- e r(i + c ) 



rn 
s 



(17) 

(18) 
(19) 



3.2 Solution obtained by the application of IBFE+IBP 

The first step required for combining these two techniques is to eliminate the mass from the diagram in (j4|). 
Once we have done this, we can apply IBP to the massless triangle @. For this purpose, we will use IBFE on 
the propagator that contains this mass and expand it according to this technique, as it will be shown next: 



1 a f 2\»i i 2\«2 iy + n\ + n 2 ) 



(q 2 — to 2 ) 



or in a graphically equivalent way 



yq z — m l ) v 



u (u + ni + n 2 ) 

WT 



(~n 2 ) 



Therefore, (HJ will be given by 




<3 = ^"i. 



2\m ( a i n 2 ) 

-to —, — x 



r(oi 




(20) 



(21) 



(22) 



The next step is obtained from the application of IBP to the triangle, once the mass has been eliminated from 
the graph. However, before this, and considering the invariants in (j4]) we can rewrite the triangle rule ([2]) as 




In order to apply this rule, we will consider special values for the indices: we take 03 = 1 and keep arbitrary 
values for ai and a 2 (04 = a 5 = 0, since they are external lines). Then, the previous expression can be written 
as: 



G 



Applying (gU) in we have: 



m ,a. 




E 0ni,n a (~m 2 ) 



E 



(-n 2 ) 




(25) 



For simplicity, let us define: 



-n 9 + 1 



A =-H 



. m (oi + ni + n 2 ) T(D + n 2 ~ a 2 - 2) T(n 2 + 1) 
' f>0 r(£ + n 2 - a 2 - 1) r(n 2 ) 




(26) 



and also: 



T(a 2 + 1) 



r(a 2 ) ^ Vni '" 2 



/ 2 \ ™i ( a i + n i + n z) r(-D + n 2 - a 2 - 2) 
' T(ai) r(£> + ra 2 - a 2 - 1) X 




(27) 



where the identity tr 



r(a + i) 
r(a) 



has been used to express everything in terms of Gamma functions. 



To find the MRE's associated to the terms A and B, we only need to express the bubble diagrams in these 
equations through the following identity [T5] : 
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a 




p =J2GA(a,P; nj )x { ^ p , (28) 



where the one bubble function is given by 



n ( r \- t- 1 )"^ V" ^ A!A 2 A 3 

Ga ^ " r(a)r(/3) 2- ^,„ 0+1) ,n W+a) ,d y ( 29 ) 



the constraints being {A^} 



A l = ( 2 + n J + + n U+2)) , 

A 2 = (a + n,j + , (30) 

A 3 = (j8 + rij + ny +2 )) ■ 



Through a careful replacement in A and -B we find for A 



A 



(-I)" 



r(£> + n 2 - a 2 - 2)r(ra 2 + 1) 



r(oi)r(a 2 ) ni ^ TlB , "* 1 '-'* B v ■" ' v ~' r(£> + n 2 -a 2 -i)r(n 2 )r(i-n 2 )r(f + n 3 ) (31) 

x (ai + m + n 2 ) ( y + n 3 + n 4 + ri 5 ) (1 — n 2 + n 3 + n 4 ) (a 2 + n 3 + n 5 ) 
and _B is given by 

(-1)"* ^ x , _.2\"i,.™ T(D + n 2 -a 2 -2) 

D , _ \ 

(32) 



B 



r( ai )r(a 2 ) 



E <t>m,..,n B (~m 2 ) 1 (s) 



r(D + n 2 - a 2 - l)r(-n 2 )r(f + n 3 ) 



x (ax + n x + n 2 ) + n 3 + n A + rig) {-n 2 + n 3 + n 4 ) (1 + a 2 + n 3 + n 5 ) . 
Once we have obtained the MRE's of the terms A and B, we proceed to find the corresponding solutions 



according to the quotient 



rn 
s 



3.2.1 Solutions from the term A 



Solution for the region 

region: 



< 1 : From the MRE pi]) we get the following solution for this kinematical 



At 



2 + a 1 +a 2 - D , 2 + a x + a 2 - D , 1 + a x + a 2 - § 
i + ai+a 2 -D , 2 + ai - f 



m 
s 



where the factor rj is given by : 

d ^ _ ai _ 02 _! r (oi + a 2 + 1 - § ) r (f - a 2 ) r (f - ai - 1) r (£) - ai - aa - 2) 



r? = (-l)-2 ( a ) 2 
and also the term 



r( ai )r(a 2 )r(^- ai -a 2 -i) z 



.4- 



r , a 2 , 1 + a 2 - f , 1 + a 2 - f 
2 + a 2 -f 2 , f- ai 



m 
s 



(33) 



(34) 



(35) 
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where we have defined r? as 



Solution for the region 

A 3 

where 



> 1 : In this region we only have one term : 

d 2 • ! r (I ■ (l-> ~r -r - 1 



= V x 3-F2 



1 + ai + a 2 - ■ 

2 ' 2 



riv 



(37) 



?7 = -(-1) 2 (-TO 2 ) 



I) 



-1 r (1 + ai + a 2 - f ) r (f - a 2 - 1) r (f - 1) r (2 + a 2 - f ) 
r( ai )r(i + a 2 -f)r(f) 2 



(38) 



3.2.2 Solutions from the term B 



Solution for the region 



< 1 : For B we have as solution in this kinematical region 



/ 2\ / 2 + a 1 +a 2 -D , 2 + a x + a 2 - D , 1 + ai + a 2 - 

M^J =??>< ^{ 3 + a 1+ a 2 -D , 1 + a, - § 



TO 



being 



ri = (-l)"f ( s )T- a i- a *- 



! T (ai + 02 + 1 - f ) r (f - a 2 - 1) T (f - 01) T (£) - ai - a 2 - 2) 



and also the term 



5 2 = r? x 3 F 2 ( 



r(o 1 )r(o 2 )r(£)-ai-a2-l) 



1 + a 2 , 2 + a 2 - % , 2 + a 2 1 



.5 + a 2 - T , 3- 



2 _ ' ^ 2 

1 — ai 



m 
s 



where 



(39) 



(40) 



(41) 



^=(-l)-t( a )- 1 -«" (-TO 2 ) 2 



§ - Bl r (ai - § ) r (f - 2 - a 2 ) r (a 2 + 1) 



r(oi)r(a 2 )r(f -i-a 2 ) 



(42) 



Solution for the region 



> 1 : In this case we only get one hypergeometric function: 



B 3 (^)= V x 3 F 2 ( 



1 + a 2 , 1 + 01 + a 2 - # , £ - 1 



2 ' 2 



2 ' 2 



TO 



where 



(43) 



n = (-!)" 2 (-m 2 ) 



D 



1 r (1 + ai + a 2 - f ) r (f - a 2 - 1) r (f - 1) r (1 + a 2 ) 
r( fll )r(a 2 )r(f) 2 



(44) 
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3.3 Comparison of the solutions for the region 



77T 



< 1 (unitary indices) 



In order to compare the solutions we have obtained with both techniques, we will take the result for the 
* <1. 



kincmatical region 

For the case of unitary indices and D = 4 — 2e, the combined technique IBP+IBFE produces 

G(^)=B 3 (^)+A 3 (^), (45) 

i.e. 



G 



IBP+IBFE 



te)= -(-l) e (-™ 2 ) 



_i_ £ r(l + e) 
e(l-e) 2 



3F2 



2 , 1 + e , 1-6 
2-e , 2-e 



(46) 



-e 3F2 



1 , 1+e , 1-e 

2-e , 2-e 



In the same region of interest, the solution given in flo]) obtained by a unique application of IBFE is notoriously 
different, both from the amount of terms as well as from the structure. Let us discuss this point more in detail: 



G 



IBFE 



(^)=-(-l) £ (-- 2 ) 



-1-e r(l + e) 
e(l-e) 



1*1 



1 + e , 1 

2-e 



(47) 



We can do now comparisons to show the equivalence between both techniques (see the Appendix for a general 
proof). A direct possibility is to analyze the equivalence in certain limit cases, i.e. to compare, for m = the 



solutions we obtained from the expressions valid in the kinematical region 



< 1 and to compare also the 



solutions for s = 0, this time considering the corresponding expressions for the region 

3.3.1 Analytic equivalence in the limits m — > and s — > 
Solution for m = : In this limit, the Eq. (|46|) has the form 



rn 
s 



> 1. 



GiBP+iBF E \ m=0 = A 1 (0) + Bi (0) = 2(-l) e ( s y 



r(i + e)r(i-er 

2e 2 r(l-2e) 



whereas ([4"7]) 



_ 1 - e r(i + e )r(i-e) J 



e 2 r(l -2e) 



GlBFE\ m=0 = (-lY(sy 

Some algebra allows us to show that effectively, 

GlBP+IBFE I m =0 _ G IB FE\ m=Q = 0. 

Solution for s = : The Eq. (|4"o]l becomes now 

Gibp+ibfe\ s=0 = A 3 (0) + B 3 (0) - - (-l) e (-m 2 )^ 1 ^ ^7^4- 
In the same way as (1471) . in this limit we have: 



^ i , ,,«/ 2\-!- e r (1 + e) 

GlBFE\ s=0 = -(-iy(- m ) ~ (1 _ e) ■ 



(48) 



(49) 



(50) 



(51) 



(52) 
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In this case both solutions are equivalent: 

Gibp+ibfe\ s= o — Gibfe\ s=0 = 0. (53) 

According to the results we have obtained previously, we can say that IBP+IBFE is a technique that can be 
used in an equivalent way to IBFE to evaluate Feynman diagrams. Nevertheless, the former has some technical 
advantages that extend the usefulness of IBFE beyond the family of diagrams, where it can ideally be applied 
and which is described in [7]. 



4 A two loop application 

As an example in a more complex scenario, we will use now the combined method IBP+IBFE for computing 
the correction to the propagator given by two loops and five propagators. It is possible to evaluate such a diagram 
for the massless case and with integer indices in their propagators iterating the use of IBP (if it is necessary) , 
transforming the original topology into a sum of topologies which, for the unitary indices case, generates more 
simple graphs which can be reduced loop by loop in terms of bubble diagrams. In this work, we will evaluate 
this diagram taking now one massive propagator, mass M [14] . with index a\, as we can see in the following 
graphic equation: 




region of interest. We impose the following conditions to evaluate the diagram: 



• The indices 03,04, and 05 are integer quantities. In particular, we will have 03 = 04 = 05 = 1. This 
assumption does not imply a loss of generality, since the indices in the loop integrals are unitary. 

• We will consider arbitrary indices a\ and a-i- This will allow us not only the possibility of taking several 
values for the indices, but also to implement topological variations of the diagram by inserting corrections 
to the propagators, having then more loops, or by taking equal or different masses for these propagators. 
The last point is related to the fact that the extraction of the mass from a certain propagator generates a 
massless propagator with a variable index, as it is shown in the Eq. (|2ip . 

• Let us assume the masses 7712 = 7713 = TO4 = — and m\ = M. 



Let us then consider the loop integral for this graph: 

d D qi d D q 2 1 



G 



1 



1 



i7T D/2 ijT D/2 ( g 2 _ M 2) 0l 



(qi+pY 



(qa + Qi) I2 (p-qzY 



(55) 



The idea is to extract the mass M from the integral in such a way that the resulting integral corresponds to a 
massless diagram and then, in a second step, to apply the triangle rule ©. 
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This is the place where IBFE can be used. This technique allows us to find the MRE of the propagator that 
contains the mass , i.e. 



1 _ V 6 (a 2 Y 11 (-M 2 ) n2 + n i +"2) 



(56) 



By replacing then in the integral (|55p . we get the following expansion (In this expansion, we must keep in mind 
that the sum indices are not necessarily integers): 



r _ a i , f2 \"2 (ai +^1 +"2) f d D q 1 d D q 2 1 1_ 

it, [ ' r ^ J i- D/2 ^ D/2 (ii)- ni [( qi+P y 

We can represent also in a graphic way this result in the following way: 



111 



(92 + 9i) 2 <?2 (p-fla) 



\2 ' 



(57) 




At this point, we can apply the triangle rule to the resulting topology. We have to mention that the multiregion 
expansion of the massive propagator contains in its structure the expansion for the cases of small and big masses 
simultaneously, which is precisely what we need to compute this diagram in all of the interesting kinematical 
zones. 



4.1 Application of the triangle rule 

To apply the formula we will consider the left triangle of the diagram in (|58p. As a result of this, we get in 
a simple way the following expression for G: 
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G= E 



-M 



2x112 (oi +ni+ n 2 ) 



r(ai) (U + m-oa-2) 




(59) 

We will express all the coefficients (only those which contain some sum index) of the previous equation in terms 
of Gamma functions. This will be useful to form the resulting hypergeometric functions. We use the following 
identity: 



r (a + 1) = aT (a) => a = 



r(a + l) 



r» ' 

with this, we may replace in l)59p the sum index dependent coefficients 

r(-m + i) 



(-m) = 



and for the denominator, 



1 



r(-m) ' 

r (D + ?n -o 2 -2) 



(D + ni - a 2 - 2) T (D + n x - o 2 - 1) 

4.2 Analytical Solution 

According to the Eq. (|59"j) . we can write the original diagram as a sum of four terms, i.e. 



(60) 



(61) 



(62) 



G = At + A 2 + A 3 + A 4 



(63) 



which are defined as 
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M= E ni , na (-M 



. , ,,,2 (qi + ni + n 2 ) T(D + ni ~ a 2 -2) T {-ni + 1) 

r(ai) r(D + m-a 2 -i) r(-m) 



•"7+1 




(64) 



A 2 = E &»i,n a (-M 2 ) 



2 , „ 2 (ai + ni + n 2 > r (D + m - a 2 - 2) T (ni + 1) 



r(oi) r(D + m-a 2 -i) r(m) 



■"7+1 




^3= (02) E 



-M 



2 \n 2 (ai + ni + n 2 ) T (£> + ni - a 2 - 2) 
r(ai) r(L> + n 1 -a 2 -l) 



(65) 



and, finally, 




a 2 + 1 



(66) 



M = (-02) E 0m,n 2 
m ,712 



(-M 2 ) 



2\n 2 (ai + ni + n 2 ) r (D + m - a 2 - 2) 



r(ai) r(D + m-02-i) 




(67) 



a 2 + 1 

The resulting topologies in each term Ai (i = 1, 2, 3, 4) are appropriate to compute them directly with the IBFE 
integration technique. With this technique, we will be able to get the analytical solution in the possible kinemat- 
ical regimes where each term Ai contributes. In what follows, it is necessary to get the equivalent complete MRE 
for each one of these terms to obtain then the analytical solution. We have two ways to get this expansion: the 
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first one is to expand the parametrical integral that corresponds to the complete diagram present in each term 
Ai, and the second one, which is the most straightforward procedure, is to use the idea of modular reduction of 
the diagram [TS] , that allows us to find very rapidly the MRE of each one of the diagrams as a product of MRE 
of their constituent individual loops. These are known as loop functions. We will use this last form to start with 
the analysis of the remnant diagrams in each term Ai. We will present the complete calculation procedure for 
the first term of (|59p . A\, and then the other terms will be presented in a summarized way. 



4.2.1 Complete solution of Ax 

To find the solution from this term, we need first the MRE associated to the diagram. We will use the idea of 
loop function representations [15] to obtain the expansion which has the form : 



■n 1+ 1 




P = y"]GU(-ni + l,a 2 ; ra 3 ) ^<3a(1, 1; n 6 



x (-"3-"«) ,( 68 ) 
• p 



where the loop function Ga corresponds to the MRE of each loop, or part of the diagram, and which has the 
following MRE: 



Ga (a j7 a k ; n,) 



(-l)-s- AiA 2 A 3 

r(a J )r(a fc )„,. i ^_^' f,( * flJ ' Wt ^ m + mY 



(69) 



«(i+l),™(i + 2) 

In this case the constraints {A^} are given by the following expressions: 

A l = ( f + n, + n (l+1) + n (t+2) ) , 
A 2 = Uj + rii + n {l+l) ) , 
A 3 = (cik +rii + ri( l+2 )) • 



(70) 



Clearly, we have to replace the pictorial representation of the propagator by its corresponding algebraic expres- 
sion: 



1 



P (pa)"' 



Some algebra allows us to find the MRE of the diagram in ((68 



(71) 



-n 1+ 1 




(-1)- 



r(-m + i)r(a 2 ) 



(72) 



x E 



[p 



2\™3+™6 



A (l) A (l) 



r(f + n 3 )r(f + n 6 y 



whose constraints are determined as follows: 
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f rc.4 H 


-"5) , 


A« 


= <-ni + 1 


+ TI3 - 


f "4) , 


A« 


= (»2 + n 3 - 


f n 5 ) 




A« 


= (# + n 6 - 


f n 7 4 


- «s) , 


A« 


= (l + n 6 + 


"7) , 




A« 


= (l + n 6 + 


«s) ■ 





(73) 



Finally, by introducing the necessarily replacements, we have obtained the MRE for the first term of the solution 
of the diagram G. We have added a superindex in the constraints. This only refers to the fact that these 
constraints are associated to the first term of G, A\. So, we get then the MRE for this term: 



A 1 



(-1)- 



r(oi)r(o 2 ) 



-M 



2\™2 / 2\™3+«6 



T(D + ni - a 2 - 2) 



Af»...A«A 



r (D + ni - a 2 - 1) r (-m) r(f + n 3 )r(f + n 6 ) 



where we have also defined an extra constraint associated to the MRE of the massive propagator: 



(74) 



A = (ai + ni + n 2 ) 



(75) 



Starting from the expression (f74|) . it is possible to find the solutions in the two possible kincmatical regions 



M 2 



< 1 and 



M 2 



< 1, where the solutions for each kinematical region can be divided into groups, as a sum 

of generalized hypergeometric functions and where the sets will be defined by the terms A\ ( jj^j and A\ {^r^j > 
respectively. 

We will show now the analytical solutions we found for the term A\ according to the corresponding kinematical 



regions: 



Solution in the region 



M 2 



< 1 



, /V \ v p ( 1 + ai + a 2 - f , 
A i \ m* ) = V x 3-F2 ^ d 



D 

2 



1 



P 

M 2 



where the factor 77 is defined as: 

D -2 f °-i- ai - a2 r (1 + ox + « 2 - f ) r (f - 02) r (2 - g ) r (f - 1) 



r, = (-l)- D (p*) 2 -*(-M 2 y 



r( 0l )rp-2)r(f) 



(76) 



(77) 



Solution in the region 



M 2 



< 1 



Ar(f)= ,X 3 F 2 ( 2 + a 1+a2 -D fl ,_2 + fll -» 2 - Z> 



1 + ai + a 2 - y 
, 3 + ai + a 2 — D 



+ V x 3 F 2 



1 + 02 - f , 1 + o 2 - f 
§ — ai , 2 + a 2 — 4r 



«2 



M 2 
p 2 



M 2 

p2 



(78) 



being the factors 
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r (l + ai + a 2 



D 



D 



, r , 



a 2 )r( 2 -f)r( 



' D 
2 



ai-1 



r(oi)r(o2)r(D-2) 



(79) 



rp- ai -a 2 -2)r(f -l) 



and 



R = (-1)-° (p 2 ) 



T{D-ax-a 2 -l) 

d 1 _ air(1 + ai _ f) r(f- a2 -i)r(2-f)r(f-i) 1 



(80) 



r( ai )r(D-2)r(f -a 2 ) 

Nevertheless, the solution of the diagram G is by no means complete, since we have to consider also the terms 
that can be generated from the MRE associated to {^2,^3,^4}. Then, we have to add them algebraically, 
according to the kinematical region of interest, i.e.: 



G 



A, (£) + A 2 (£) + A 3 (£) + A, (£) 



for M 2 > \p 2 



A x 



A, 



A. 



3 (^)+A 4 (f ) , for \p 2 \>M 2 



(81) 



We will proceed in an analogous way, as we did for the calculation of Ax, for the other terms A ; : 

4.2.2 Contributions to the solution generated by A 2 

For the diagram that appears in A 2 , the MRE is obtained as 1: 



■n 1+ 1 



_P = y~] Ga(1, 1; n 3 ) ^2 G A {-ni + l,a 2 - n 3 ; n 6 ) x 

"3 ™6 



-n 6 



, (82) 



and we use then the expression (|69[) to replace the one-loop functions G a in the expansion of this graph, getting 
in an explicit way the MRE of the term A 2 : 



Ao 



r( a l) m,...,n. 



M 2 )" 2 (p 2 )" 6 r — + 711 - 0,2 - ^ r ^ + — 



r(U + ni-oa-l)r(ni) 



AfL.A^A 



T(f + n 3 )r(f + n 6 )T (-nx + 1) I> 2 - n 3 ) ' 
where the constraints correspond to the following expressions: 



(83) 





= <f + n 3 - 


f n 4 4 


-"-5) , 


A 2 2 > 


= (1 + n 3 + 


nx) , 




Af 


= (l + n 3 + 


"-5) , 




Af 


= (f +n 6 - 


f n 7 4 


-«s) , 


Af 


= (-"1 + 1 


+ n 6 - 


f n 7 >, 


A (2) 
^6 


= (a 2 - n 3 - 


f n 6 4 


-n 8 ) • 



With this information, we can find in this case the solutions for both relevant kinematical regions: 



(84) 
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Solution in the region 



Mi 



M 2 



< 1 



3 + a x + a 2 - L> , 2 + a 2 - f , 1 



2 ' 



m 2 y ' 



(85) 



where the factor rj is given by the following expression: 



^ = { _ iyD ^_ M ^D-3- ai -a 2 r(£>-3- fl 2 ) r (3 + ai + a 2 - £)) T (4 + a 2 - D) V (2 - § ) T (§ - l)' 



r (01) r (£) - 2) r (f ) r (3 + a 2 - d) 



(86) 



Solution in the region 



M 2 



< 1 



{£)= ^ 

+ 77 x 2 Fi 



4 + ai + a 2 - ^ , 2 + ai + a 2 - L> 



2 + ai - f 
3 + a 2 -£> , 1 + a 2 - § 



m 2 

p2 



M 2 

p2 



where the factors are defined as 

.D-3-ax-a, r (3 + oi + a 2 - £)) T {D - a 2 - 2) T (2 - § ) T ( § - 01 - l) 



n= (-i)- D (p 2 y 



r(ai + l)r(-ai)r(£)-2) 



r (D - ai - a 2 - 2) T (§ - l) 2 r (1 - < 



' T (D - ai - 02 - 1) r - ai - as - 3) r (2 + 02 - f ) ' 



I) 



r7 = (-l)- c (p 2 )2— (-M 2 ) 2 



» 1 _ ai r(i + o 1 -f)r(f -a 2 _i)r(2-f) 2 r(f-i) £ 



r( ai )r(/)-2)r(f-a 2 )r(f)r(i-f) 



(87) 



(88) 



(89) 



4.2.3 Contributions to the solution generated by A 3 

We will follow the same procedure for getting this term, i.e. to find the MRE of A 3 and their corresponding 
analytical solutions: 




V = y^GU(-ni,o 2 + l; rc 3 )^GU(l,l; n 6 ) 



(-n 3 - n 6 ) 



(90) 



a 2 + 1 



(— M 2 ")™ 2 ( 2N t™ 3+ ™ 6 ^ (-D + ni a 2 2) ^ ...^-» 6 

1 ' [p ' r (£> + m - a 2 - i) r(f + n 3 )r(f + n 6 )r (-m) ' 



A< 3) ...AfA 



r(ai)r(a 2 ) 



with the constraints: 



n 1 ,...,n 8 



(91) 
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A? 
A 2 3 

A? 



= (f +n 3 + n 4 + n 5 ) . 
= (-m + n 3 + n 4 ) , 
= (a 2 + 1 + n 3 + n 5 ) , 
= (f +n 6 + n 7 + n 8 ) . 
= (1 + n 6 + n 7 ) , 



Ar = (l + n e +n s ). 



Solution in the region 



< 1 



M 2 

A 3 (%f)=VX 3 F 2 



1 + ai + a 2 - 



2 

D 

2 



f - 1 , 1 + a 2 

D 
2 



M 2 



77 = (-1)--° (-M 2 ) 2 
T (l + ai + a 2 



-1 — Ol— 02 / 2 



D 



^r(f-a 2 -i)r(i + a 2 )r(2-f)r(f-ir 



r(ai)r(a2)r(D-2)r(§; 



Solution in the region 



M 2 



< 1 



, , 2 + ai +a 2 -D , 2 + oi + a 2 - D , 1 + Oi + a 2 - f 
= n x 3-F 2 -j | Do,, n 



M 2 
p 2 



1 — ai + y ! o + a 2 — 



2 + a 2 - f , 1 + a 2 



M 2 " 

P 2 ' 



where the different factors are given by 
V = {-l)- D (p 



D . 2 ^_ 3 - 01 -a 2 r (£> - ai - a 2 - 2)r (f - o 2 - 1) r (2 - f ) r (f - 1)' 



r(ai)r(a 2 )r(£)-2) 



x r(f - ai )r(i + a 1 + a2 -f) 

T(D- ai -a 2 - l) 2 

„ = (-1)"" (p 2 ) ^ 3 - a2 (-M 2 ) ^ #) T (1 + a 2 ) r (2 - f ) T (f - a 2 - 2) T (f - l) 5 



r(a 1 )r(ai)r( J D-2)r f -o 2 -i 



4.2.4 Contributions to the solution generated by A4 



P = *^2 G ' 4 ( 1 ' 1 ! n 3) X! _ n 3; fl 2 + 1; «6) X 



-n 6 
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(-1) 



-D 



r(ai)r(o 2 ) 



E 



(-M 2 )™ 2 (p 2 )" 6 r + ni _ ° 2 _ 2 ) 



A^...Af A 



r (d + m - a 2 - 1) r(f + n 3 )r(f + n 6 )r (-m - n 3 ) : 



Al 4 

Ar 

Ai 4 
A^ 4 
A^ 4 



(§ +n 3 +n 4 + n a ) , 

(1 + n 3 + n 4 ) , 

(1 + n 3 + n 5 ) , 

(f + «e + n 7 + n$) , 

(—m -n 3 + n e + n 7 ) , 

(a 2 + 1 + n e + n s ) . 



Solution in the region 



p- 



M 2 

A 4 (fe) = V X 3^2 



< 1 



D — 2 , ai 
D - a 2 - 1 



3D 
' 2 



£-1 

a 2 - 3 



P 

M 2 



. r , 2 + a 2 - f , 3 + a 2 + fll - D , a 2 + 1 
+ ??X 4F3( f , 2 , 4-D + a 2 



M 2 J ' 



being 



v = (-iy D ( P 2 ) D - 3 - a2 (-M 2 y ai 



r (3 + a 2 - -D) r (f - a 2 - 1) r (£> - a 2 - 2) r (f - 1) 

f M r(f-d 2 -3)r(fl- fl2 -i) 



-D / 7i^2\- D -3-ai-a2 



5 =(-l)-»(-M 2 ) 
Solution in the region 



r (£> - a 2 - 3) r (3 + a 2 + ai - £>) r (2 - f ) r (f - 1) r (1 + a 2 ) 

r( ai )r(a 2 )r(D-2)r(f) 



M 2 



< 1 



r? x 3 F 2 
+ V x 3 F 2 



4 + ai + a 2 - ^ , 2 + a 1 +a 2 - D , a x 



2 + ai-f , 3 + ai-Z) 



M 2 

p2 



2 + a 2 - f , D - 2 , a 2 
D - ai - 1 , § 



M 2 

^2 



with 



/ n -D / 2\P-3-»i-» 3 r (f - a 2 - 1) T (£> - ai - 2) T (£> - ai - a 2 - 2) 
71 ~ ( ' ^ > " ' T(a 2 )T(D-2) 

r (2 - f ) r (3 + ax + a 2 - D) V (| - l) 2 
X T (D - oi - a 2 - 1) T - ax - a 2 - 3) T (2 + oi - § ) ' 

^=(-l)^( P 2 )- 1 - a2 (-M 2 ) 



D/2X-1-2/ ll ^o- 1 -2r(2 + ai-D)r(2-f)r(f -l) r(i + « 2 )r(-a 2 ) 



r( ai )r(a 2 )r(f)r(i-a 2 ) 



(99) 



(100) 



(101) 



(102) 



(103) 



(104) 



(105) 



(106) 
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4.3 A particular case: Solution for M = with indices a\ 



a 2 



In what follows, we will test the solution we got for the diagram (|54[) . For this, we will compute the massless 
case, which has a trivial solution when integration by parts is used. In this case, we have to consider the solution 



obtained by IBP+IBFE in the region 
G(M=o) = A 1 



P 



< 1 and to take it's limit M = 0, i.e. 



.4. 



A. 



A 4 



M=0 



(107) 



It is easy to verify that in this limit only four terms will survive. Using some algebra, we can get the following 
result: 



G 



IBFE, 



2 (-1)"^ (p 2 ) 



D-5 



r(f-2)r(f -i) 2 r(2-f)r(i?-4) 



r(f-i)r( 3 -f) 

T(D- 3) 2 



T(D-2) 
T(5-D) 



r(f -s)r(3-f) 



(108) 



On the other side, the conditions under which we now have to evaluate the diagram are such that the solution 
can be found in a conventional way, i.e. applying IBP through the triangle rule. So, we find 



(109) 




{D-A) 





and we can rewrite this as 



G 



IBP,, 



1 G(l,l)[G(2,l)-G(2,3-§)] ( P y- 5 . 



(M=0) (£> - 4) 
where the factors G(a, /3) are given by 

Finally, we obtain the following result for this two-loop diagram: 



(110) 



(111) 



Gibp (m=0) = (-1)- D (P 2 ) 



D-5 



D 



i r 



' D 



{D-A) 



r(D-2) 



i)r(3-f) r(5-D)r(£)-3) 



r(L>-3) 



r(f -5)r( 3 -f 



(112) 



The solutions (|108p and (|112p seem to be apparently different, but by using the properties of the Gamma function 
and some algebra, it is possible to show that both terms are equal, i.e. 



G 



IBP,,, r=n 



G 



~(M=o) " IBFE( M=Q } 

This result gives support to the most general solution for this diagram in the case of one massive propagator 



0. 



(113) 
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5 Conclusions and comments 



It is not always possible to combine the IBFE technique with another conventional integration method. For 
example, let us evaluate the one-loop diagram where one of the propagators has a mass M: 



G = 



d D k 



1 



(114) 



in D/2 (fc2_M 2 )°i [(p-fc) 2 ] Q2 ' 

As a first step, we will take out the mass M of the massive propagator in such a way that we may associate the 
resulting integral to a massless diagram. We need to express this propagator in terms of a MRE, i.e. 



W = E ( k T (-m 



2 \n 2 (oi + ni + n 2 ) 



(k 2 - M 2 ) 

finding, when replacing this result in (|114[) . 

2 \n 2 (ai + ni + n 2 ) f d D k 



r(oi: 



G 



y ] 4>n 1 ,n 2 
ni,n 2 



-M 



r( 0l ) 



i?r D/2 (fe2)-m [(p- jfc)2p 



(115) 



(116) 



If the idea is to use only IBFE , the next step would be to find the corresponding MRE for the resulting loop 
integral. However, instead of doing this we will evaluate directly, and in a conventional way, the corresponding 
integral. This integral is elementary and represents a massless bubble diagram whose solution is given by 



d D k 



(-I)' 



. fi r(-ni + 02-f)r(f + m)r. 



(12) 



i7r D/2 (fc2)-m [(p _ fc)2]°2 v r(-m)r(a 2 )r(£> + m - o a ) 

We replace this result in ()116j) getting the following MRE for G: 



(p2)-ni 



(117) 



G = (-I)"* ^fe\(p 2 )-- E ^.n, (-M 2 )" 2 b 2 )" 1 (-x + * + n 2 ) r( - ni + ° 2 - f )r(f + 



r( 0l )r(a 2 ) 



T(-n 1 )T(D + n 1 -a 2 ) 



(118) 

According to the prescription of IBFE, we notice immediately that the maximum number of possible terms we 
can extract to form the solution is 
2 

2. In fact, one of these terms corresponds to an hypergeometric series which is solution for the 



1 



kinematical region 



M 2 



P 



< 1 and the other one would give us the solution for the region 



P 



M 2 



< 1. However, 



this solution is not complete since some terms are missed [21 [3]. In fact, it remains to add one term to the 

\M 2 

< 1. The terms which can be extracted (|118p are correct, but they 



corresponding solution in the region 



P 



represent only a part of the solution, and this is the problem. 

We will find now the MRE for this diagram using IBFE as the only integration method. Then, we will 
compare with the previous result shown in (I118[) . Schwinger's parametric expansion of (|114[) is 



G = 



(-1)- 



r(oi)r(a 2 ) 



ai— 1 „a 2 — 1 



exp (xiM 2 ) exp 



dx\dx2 x^ 1 x 



X1X2 9 

T P 

Xi + x 2 



(xi + x 2 y 



(119) 



After some algebra, the following MRE for the diagram G is obtained: 

(-l)-f x . x (p 2 ) ni (-Af 2 )" 2 



G 



r(ai)r(a 2 ) 



y ^ ^ > v - ; TT a, 



(120) 
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where 



Ai = (§ + n 1 + n 3 + n 4 ) , 

A 2 = (01 + n x + n 2 + n 3 ) , (121) 
A 3 = (a 2 + ni + rt 4 ) . 

We observe that this time, the possible number of series representations that can be extracted from (|120p is 

( 4 \ 

given by the combinatory I ^ 1=4. This shows that when conventional techniques, combined with IBFE, 

are used, the MRE for a certain diagram is more simple than the MRE obtained only by means of IBFE. This 
means that the conventional techniques do not generate sums or Kroneker deltas and therefore the number of 
different ways for summing, using the constraints of the MRE, becomes also smaller. 

In general, the solutions that are obtained from these MRE's do not correspond necessarily to the correct 
solution. They can be incomplete. In summary, once we have chosen to use IBFE as an integration technique, 
the whole resulting process must be done with this technique. However, to simplify the problem, it is possible 
to apply other techniques before IBFE, but not viceversa. 

The exception to this rule is IBP, where the integration in terms of the topological cancellation of one of the 
graph lines do not produce the disappearance of terms in the solution. The explanation for this relies in the fact 
that the IBP expresses the original diagram as a sum of diagrams. In the example given at the beginning of this 
work, we used IBP in the intermediate step of getting the MRE. The result was a sum of MRE's in which each 
one contributed to the solution, with terms associated to the corresponding interesting kinematical regions. The 
conclusion is that IBP is a technique compatible with IBFE. 

In this work, we have emphasized the role of the triangle identity and the fact that, under certain particular 
mass distributions, it is possible to evaluate a series of diagrams with the procedure described above. For 
example, 




Although these examples do not necessarily generate a hypergeometric series of one variable, the procedure 
can be directly applied to these topologies. Another option for generalizing is the insertion of bubbles in the 
propagators, with arbitrary indices, which can be massless or massive. As an example, let us consider the 
following case: 




, etc. 



We have shown that IBP+IBFE is a technique that may extend the classes of diagrams that can be evaluated 
through this simple procedure, beyond the scenario given in [7] and, although here we have not considered the 
IBP in his most general version, we assume that the combination IBP+IBFE can be also useful beyond the 
triangle identity. 
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Appendix 

In what follows, we will prove that the solutions that we obtained for the one-loop example given in Eq. 
((31) : the first one using only IBFE (|47p and the second one through the combined technique IBP-IBFE (|4"6")) are 
equivalent. Basically, we have to show that the functions defined below, G a and Gp, are equal: 



and 



rq-6) 

r(2-e) 
— e 3*2 



G 



3*2 



2 , 1+e , l-e 
2-e , 2-e 



(122) 



1 , 1+e , l-e 

2-e , 2-e 



2*1 



1 + e , 1 

2 - e 



(123) 



The variable x represents the quotient ( — =• V In order to prove the equivalence of both solutions we need to 
remind the series expression for an hypergeometric function: 



9*9-1 



{«} 
{b} 



Y> («i)fc-K) fc % k 
2-,( bl)k ...( bq _ l)k k v 



(124) 



then we can rewrite (I122[) as a series, as follows 
T(l-e) 



Ga - r (2 - e) 

= T(l-e) 
T(2-e) 

On the other side, we have that 



f (2) fc (1 + e) fc (1 - e) fc x fc °° (l) fc (1 + e) fc (1 - e) k x k 
to (2-e), (2-e), fc! ' - 



(l + e) fc (l-e) fc 



k^o (2 - e), (2 - e), 
„fe" 



fc! 



(125) 



a; 



( 2 )^ r ^j fc) = (i + fe)r(i + fe) = (i + fe)(i),, 



which implies that (11251) reduces to the following equation: 



T(l-e) 



E 



q) fc (i+o fc (i-^) fc 

- (2-e), (2-e), 



[1 + fe - e] 



T(2-e) 

Then, we can transform the expression (1 + k — e) as follows: 

(1 1 fc el _ T(2 + fc-e) _ r(2-e)(2-e) fc 
( £j "r(l + fe-e)-r(l-e)(l-e) fc - 

By replacing this result in (|127[) some algebra shows finally that 



fc! 



(126) 



(127) 



Go 



(l),(l + e) fe z* 



2*1 



1 + e , 1 
2-e 



Gt 



(128) 



Although the equivalence has been shown for the kinematical region 

> 1, since both regions are related through an analytical continuation. 



< 1, it is also valid for the region 
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